Optical conductivity of bilayer graphene with and without an asymmetry gap 
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When a bilayer ol graphene is placed in a suitably configured field effect device, an asymmetry gap 
can be generated and the carrier concentration made different in each layer. This provides a tunable 
semiconducting gap, and the valence and the conductance band no longer meet at the two Dirac 
points of the graphene Brillouin zone. We calculate the optical conductivity of such a semiconductor 
with particular emphasis on the optical spectral weight redistribution brought about by changes in 
gap and chemical potential due to charging. We derive an algebraic formula for arbitrary value of 
the chemical potential for the case of the bilayer conductivity without a gap. 
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I. INTRODUCTION 

Graphene is known to exhibit special properties related 
to the Dirac nature of its quasiparticle dynamics. As an 
example, a half-integer quantum Hall effect is observedi^ 
as was predicted^i^. Bilayer graphene also possesses re- 
markable properties. When placed in a suitably config- 
ured field effect device, a tunable semiconducting gap can 
be generated with the result that the valence and the con- 
duction band no longer meet at the two Dirac points in 
the graphene Brillouin zone^i^iiiSii^. For a review of other 
remarkable properties of such systems as well as a dis- 
cussion of possible technological applications, the reader 
is referred to RefJ^iii. The optical conductivity of a 
few-layer epitaxial graphit o^^'^^ and oriented pyrolytic 
graphitei^ii^ in finite external magnetic field has been 
reported recently, as well as for graphenoi^. There have 
also been theoretical studies^iiiii^i^i^S of the conductiv- 
ity, including discussions of optical sum rule o^^'^^ which 
continue to provide useful information^^ on the electron 
dynamics. In Ref.^^ it was found that the opening of the 
asymmetry gap in bilayer graphene leads to very small 
changes in the overall optical sum. Here we consider 
the optical conductivity of such a system with particu- 
lar emphasis on the optical spectral weight redistribution 
brought about by changes in the chemical potential, due 
to charging, and to the opening of a semiconducting gap. 
In the configuration envisioned here, donor atoms are 
seeded on the upper face of a bilayer which is also placed 
in a field effect device so that the carrier imbalance in 
each layer is different as is the electrostatic potential. 

The structure of our paper is organized as follows. In 
section II, we outline the theoretical derivation of the 
optical conductivity starting from the simplest nearest- 
neighbor tight-binding Hamiltonian which includes terms 
associated with the biased bilayer configuration. In sec- 
tion HI, we present results for the case without the 
anisotropy gap, pertinent to the unbiased graphene bi- 
layer. Here, we have been able to derive a surprisingly 
simple analytical formula for the frequency-dependent 



conductivity for arbitrary choice of chemical potential 
and we have tested it against the full numerical solution 
and find good agreement, making this a very useful for- 
mula for experimentalists. In pure monolayer graphene, 
the conductivity is flat for the half-filled case of /i = 
and has a universal value of 7re^/2/i. With finite fi, the 
spectral weight below a frequency = 2fi is transferred 
to a Drude metallic response at zero frequencyi^*^ We find 
in the case of the pure bilayer, for the half-filled ^ — 
case, the conductivity retains a universal value at w = 
and at high frequency, which is now twice the monolayer 
case, i.e. ne'^/h. Structure in the conductivity is found 
at D, — J and 27, with the structure at 7 being par- 
ticularly pronounced due to a doubly degeneratre set of 
interband transitions. With finite fi, the doubly degner- 
ate set of transitions is split into two separate features 
and a new strong absorption at 7 is seen due to the nest- 
ing of two energy bands. At low energy, spectral weight 
is once again removed below 2^ and a Drude peak oc- 
curs. In section IV, we consider the case of the bilayer 
with an asymmetry gap. The main results in this case is 
the appearance of a semiconducting gap in the conduc- 
tivity for fj, = and a shift in the structure of the doubly 
degenerate interband transitions to higher energy. For 
finite /i, the metallic behavior is restored and a Drude 
peak occurs accompanied by loss of spectral weight be- 
low 2fi and the structure due to the doubly degenerate 
transitions is split, as before. The nesting feature at 7 in 
the unbiased bilayer is now shifted to a new value and is 
broadened due to imperfect nesting of the energy bands. 
Finally, we end our analysis with a discussion in Section 
V. 



II. THEORETICAL BACKGROUND 

In order to derive and discuss the optical conductivity 
of biased bilayer graphene, we need to examine the form 
of the band structure and provide an expression for the 
electronic Green's function. To this end, we begin with 
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the Hamiltoniaii for the system under consideration and 
foUow the notation given by McCann^ and Benfatto et 



alM for continuity. The single spin Hamiltonian is given 
as: 



H = ~t^{a\^^gbin + h.c.) - t^{al^b2r,+s' + h.c.) + ■y^{al^bin + h.c.) 



n.S' 
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The first two terms are the nearest neighbor hopping 
terms for electrons to move in each of the graphene 
planes, separately. The two planes are indexed by 1 
and 2, with a single graphene sheet having two inequiv- 
alent atoms labelled A and B, as the arrangement of 
carbon atoms on the two-dimensional honeycomb lattice 
provides for two atoms per unit cell. The operator bin 
annihilates an electron on the JB-atom site with site label 
n in the graphene sheet with label 1 and a\^_^_g creates 
an electron on the neighboring A-atom site in the same 
sheet positioned at n + <5, where S can be one of three 
vectors which point to the three possible nearest neigh- 
bors. These vectors are enumerated as 5 1 = — (ai-|-a2)/3, 
^2 = (2ai — a2)/3, and ^3 = — (ai — 2a2)/3, where 
ai — {a\/3/2,a/2) and a2 = (a-\/3/2, — a/2) are the unit 
vectors of the triangular sublattice for the A or B atoms, 
and a — |ai| — |a2| — VSacc with acc the distance 
between two nearest carbon atoms. Now in layering 
graphene sheets, there are several choices for stacking. 
The one under consideration here and in experiment, is 
that of the Bernal-type stacking where if the atoms are 
labelled AI and Bl in sheet 1 and A2 and B2 in sheet 
2, then the A2 atoms are stacked directly over the Bl 
atoms, but the B2 atoms are stacked over the centers of 
the Al-Bl carbon rings. Pictures of this structure can 
be found in several reference^»ii and so we do not repro- 
duce this here. Note, however, that if we use the index 
n to reference the Bl atoms, then transferring this index 
to the A2 atoms directly above gives rise to indexing the 
nearest neighbor vectors in the second sheet relative to 
the A atoms and hence the vectors are orientated differ- 
ently and labelled as 6[ — — (5^. This stacking means that 
in terms of nearest neighbors associated with interlayer 
coupling, the A2 atom is the nearest neighbour of the 
Bl atom through a direct vertical bond. Thus, in the 
Hamiltonian of Eq. the third term shows the hop- 
ping term for an electron on the Bl site to hop to the 
A2 site in plane 2. The hopping parameter is given as 7 
and is typically about 0.4 eV. There is also a possibility 
to hop from Bl to B2 or Al to A2 but these hopping en- 
ergies are very small (see Refj^ and references therein). 
Hopping from Al to B2 in the Hamiltonian is another 
possibility and its energy is larger at ^ 0.315 eV,^ how- 
ever, we do not display it in the Hamiltonian as we drop 



this term to focus on the physics and results which come 
from the main interlayer coupling, the 7 term, which is 
essential to the discussion of the bilayer configuration. 
Furthermore, we are considering a biased bilayer struc- 
ture which can give rise to the very novel physics of being 
able to tune the bilayer from metallic to semiconducting 
behavior, consequently the biasing is indicated in Eq. ([T]) 
as a lowering of the energy on graphene plane 1 by an 
amount A/2 and the raising of the energy on plane 2 by 
the same amount. This results in the last two terms of 
Eq. II]) . The quantity A is referred to as the asymmetry 
gap. The Hamiltonian transforms in the standard way2^ 
to fc-space and can be written as a matrix: 
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where (f>{'k) = —t e^^'^ — —t e^*'''^ and we have 
followed the notation of McCann^ by using an eigenvector 
^ — (oik, &2k, a2k, feik)- The energy eigenvalues of this 
matrix define the band structure. As the bilayer has 
four atoms per unit cell (doubling that of single layer 
graphene), there are now four bands and these are given 
by: 

4(k) = ^ + ^ + |</)(k)p + (-i)"r. 



r = 



A2), 



(3) 



where a = 1, 2. Of these four bands, the lower energy 
ones ±ei(k) are essentially the original graphene bands 
with low energy modification and the ±62 (k) are higher 
energy bands reflecting the dimerized bond between Bl 
and A2, which has an energy scale of 7. As the main low 
energy physics occurs at the two inequivalent K and K' 
points of the graphene Brillouin zone, the function ^(k) 
can be expanded around the K point in the continuum 
approximation (i.e., the limit of small lattice constant a) 
to be |0(k)| Ki hvp^, where vp — V3ta/(2h). With this 
approximation, if A = and 7 = 0, we would recover 
the famous graphene band structure where e = ±hvp]<. 
is the form of the dispersion around the so-called Dirac 
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points, which are two-fold degenerate for the uncoupled 
bilayer. However, if A = and 7^0, the band struc- 
ture around these points is modified to be quadratic in 
k (although still linear at larger k) and the degeneracy 
is lifted such that the dimerized bands are shifted by 7, 
as can been seen later on in our first figure. The pres- 
ence of the bias energy A, produces an energy gap in 
the band structure and a "mexican hat" structure oc- 
curs in the lower energy band ei with a minimum at 
Egi = -fA/{2^/^^^ T^) at |0(k)| ^ (/.(fc p) = hvFko, 
where fco = (A/2)^(A2 + 2-f^)/{A^ + 72) and a "hat" 
maximum at i?oi = A/2 for fc = 0. This will be discussed 
further in Section IV. This unusual band structure and 
the presence of the energy scales of 7 and A, gives rise to 
very rich structure in the frequency-dependent conduc- 
tivity, as we will see. 

With this Hamiltonian, it is straightforward to obtain 
the Green's function G through G^^ = zl — H or 
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where z — iuJm with Un — i:T(2n + 1) the fermionic 
Matsubara frequency for n = 0,±1,±2, . . . and T, the 
temperature. For our calculation of the optical conduc- 
tivity, only certain elements of the Green's function enter 
our final expression and they are: Gn, G22, G33, G44, 
Gi3, G24. As G22(A) = Gii(-A), G33(A) = G44(-A), 
and G24(A,(/)) = Gi3(— A,0*), it is sufficient to show 
only three elements explicitly: 



Gn 

G44 

Gi3 



(A ~ 2z)(A2 + 472 - 4z2) - 4|0(k)|2(A + 2; 

8(z2-,2)(^2_,2) 

(A + 2z)[(A-2z)2-4|0(k)|2] 



^) 



8(^2_,2)(^2_,2) 

(2z - A)70(k) 
2{z^-el){z^-eiy 



(6) 
(7) 



The finite frequency conductivity is calculated through 
the standard procedure of using the Kubo formulcfc^. The 
real part of the conductivity is written in terms of the re- 
tarded current-current correlation function Hai3{^ + iQ^) 
as 



Im^c,/3(^^-(-^0^ 



(8) 



where a and (3 indicate the spatial components (here we 
will be interested in the longitudinal conductivity paral- 
lel to the graphene sheets, axx)- The retarded current- 
current correlation function is also referred to as the po- 
larization function which we calculate in Matsubara for- 



malism as outlined by Mahan^^ and given as 



1/T 



dre'"-" <T,J„(r)J/5(0) >, (9) 



where r is imaginary time, is the time ordering oper- 
ator and Vm is the bosonic Matsubara frequency 27rmT 
for m — 0,±1,±2, ... and temperature T. In Eq. Q 
Ilapi^ + *0+) is the analytic continuation of Eq. ^ to 
the real axis via ivm — s- + iO+. The current opera- 
tor is the sum over sites of the site-specific paramagnetic 
current operator=2i: 



(10) 



In order to evaluate this, we require the particular cur- 
rent operator corresponding to our Hamiltonian. This 
is found via a Peierls substitution on Eq. |T]), where the 
operators associated with hopping to a new site are mod- 



ified as 



bin 



exp(-f A- r)fei„. (Note 



nn+(5 



that we have taken the velocity of light c = 1.) The 
new Hamiltonian with the vector potential A is then ex- 
panded for small A to first order and the paramagnetic 
current operator is given by (n) = —dH/d{Aa{n)). 
For our Hamiltonian and the case of considering only cur- 
rents in the graphene sheets, the paramagnetic current is 
given as 



£(n) 
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J2{sX{4^b,n+s'-h.c.). (11) 
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Note that the Hermitian conjugate piece has a minus sign 
representing a depletion of the current for hopping in the 
reverse direction. Fourier transforming to fc-space and 
summing over n, we can write the total paramagnetic 
current operator Ja{k) = — e ^j^. \1/'^i)q\[', where 



(12) 



with Vak = — (ii/^) X^al^)"'^*'' The structure of the 
equation embodies the form of j — —nev where the ve- 
locity V is proportional to the gradient of the energy. One 
can see that the factors of it and the sum over compo- 
nents of 6 weighted by e*'' '' are simply a result of taking 
the gradient of (/)(k) which is the bare energy in the ab- 
sence of A and 7, i. e., Vk = V0(k)/?i. Consequently, 
the polarization function can then be written in the usual 
bubble approximation as 
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(13) 



or, favoring the spectral function representation of the 
Green's function where 



Gij{z) — 



dcu Aij [uj') 
27r z — ' 



(14) 



we can write the real part of the conductivity as 



r 



(Pk 

(2^ 



(15) 



where f{x) = l/[exp(x/r) + 1] is the Fermi function and 
/i is the chemical potential. In presenting our results, 
we evaluate this equation at T=0. Taking the trace and 
dropping those terms which vanish in the averaging over 
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momentum, the structure of the equation to be evaluated 
reduces to knowing the three spectral functions associ- 
ated with the three Green's functions mentioned earlier. 
The longitudinal conductivity axx{^) = cr(ri) becomes 



2n 

+Aii{ui,-A)Aii{uj 
+2[Ai3(c^,A)A*3(co 



f^[/(c.-Ai)-./(c. + f^-/i)] 



/cPk 
j^\vk\'^{Aii{u;, A)A44(w + n,A)+ Aii{ij, A)Aii{u; + Q, A) 



(2^)2 

A) + Aii{Lj, ~A)Aii{uj + n, -A) 



n, ~A) + Al^{uj, ~A)Ai3{lu + n, A)]} . 

I 



(16) 



We evaluate this in the continuum approximation around 
the K point of the graphene Brillouin zone, where |wkP = 
Vp and |(/'(k)| = Tivpk and the integral over k has a 
large upper cutoff typical of the large bandwidth. Thus, 
we introduce a factor Nf — 4 which comes from a sum 
over spin (not included explicitly up till now) and a sum 
over the two inequivalent K points {K and K') in the 
graphene Brillouin zone. Previous to Eq. (|16p . the in- 
tegral over <Pk was to be taken over the Brillouin zone, 
but in Eq. (fTB|) it is to be interpretted as over a single 
K point which would be a cone for the decoupled bilayer 
graphene. If we take the limit of this expression (Eq (HH])) 
for 7 = A = /i = 0, we find a constant as a function of 
frequency given as cfq = e^/2?i, which is twice the result 
for single layer graphene^'*. We will refer to this as the 
conductivity of the uncoupled graphene bilayer. Finally, 
for illustration, we write an example of the spectral func- 
tions here: 



where 



Ai3(w,A)= [ai3(a, A)(5(w-ea)-t-ai3(a,-A)(5(w+e„)], 

(17) 



Q=l,2 



oi3(a,A) = (-1) 



„7r (2e„- A)70(k) 



(18) 



The other spectral functions follow from the Green's 
functions in a similar manner. For numerical work, we 
write the delta functions in the spectral functions as 
Lorentzians with a broadening 77, i.e., 5{x) — {rj / tt) / [x^ + 
77^] . In the optical conductivity, this manifests itself as an 
effective transport scattering rate of l/rimp — 277 due to 
the convolution of two delta functions in the conductivity 
formula. 



III. RESULTS WITHOUT ANISOTROPY GAP 

At zero temperature, for the case of A = 0, it is pos- 
sible to derive a closed algebraic formula for the bilayer 
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conductivity. It has the form: 



2(1] + 7) 2(r!-7) 

2 



ein - 27) 



29? 



: [e(f] - 2^ - 7) + e(f] - 2/1 + 7)]e(i7 - 7) 

+ a{^i)5{9)+h{^i)5{Vl~-f) (19) 



with 



4Ma^ + 7)^4M,.-7)q(^_^^ 



2/i + 7 



and 



hi ^ ^ — In — — -^0(/i — 7) 



7 



7 



(20) 



(21) 



where ctq = I2h which is twice the conductivity of a 
single graphene sheet. This expression correctly reduces 
to the form given by Abergel and Farke^* in the limit of 
/i = 0. In the top frame of Fig. [2 we show results for 
CT(ri)/CTo as a function of ri/7 for the original Abergel and 
Fal'ko case of = (blue dashed curve), thus reproduc- 
ing their results, and our extension to finite [i^ specifically 
/i = 0.27 (solid red curve) and /i = I.27 (dash-dotted 
blue). In the two last cases, the two delta functions at 

= and 17 = 7 are shown as vertical arrows. Their 
weight is given by a(^)/2 and h{\i) of Eqs. (|20|) and ((2T|) . 
respectively. Note, that only half of the delta function at 

= is to be assoicated with the optical spectral weight 
for > 0. The arrows in Fig.[T^ are shown schematically 
to represent their relative weight. The quantities a(/i)/2 
and &(/Lt) are shown in Fig. [5] as the long-dashed blue 
and short-dashed red curves, respectively, as a function 
of chemical potential /i normalized to the plane-to-plane 
hopping 7. Also, shown for comparison is the case of 
the uncoupled graphene bilayer with 7 = but /i finite 
(solid black curve). We see that the amount of spectral 
weight in the Drude delta function centered at 12 = is 
less than it is for pure graphene, except in the limit of 
/i ^ 0, when they are equal. 

Besides the two delta functions, the finite /i result has 
regions of optical spectral weight lost when compared 
with the /i = case. In the solid red curve of Fig. [1] (top 
frame), the region from to 2/i is completely depleted 
while a second region of partial depletion is seen above 
17 = 7. The optical spectral weight lost below 2/u, relative 
to the /i = finite gamma case, can easily be computed 
fromEq. (HH) and is c(/x) = /i-|-(7/2) ln[(2/i-|-7)/7], which 
is shown as the dash-dotted green curve in Fig. O For 
finite /i (/i < 7), c(/i) > a(/i)/2 so that part of the lost 
spectral weight has been transferred to the second delta 
function at 17 = 7. The existence of the two regions of 
depletion just described is encoded in the theta functions 
of Eq. (fT9|) and is easily understood in terms of the en- 
ergy dispersion curve diagrams displayed in the middle 
and bottom frames of Fig.[l] where we show possible op- 
tical transitions. We start with the middle frame where 



the chemical potential, shown as the horizontal dotted 
line, is set at /i = 0. The transitions are vertical and all 
connect a filled valence band to empty conduction band 
states. From left to right, the arrows depict typical tran- 
sitions for the four possible ways of connecting the four 
bands: from —£2 ^ +ei and — ei — > -|-e2 , which are both 
restricted to photon energies 17 > 7, and —t\ — > -|-ei and 
— e2 +^2, with the latter one restricted to 17 > 27. 
The onsets of transitions at 17 = 7 and 27 are clearly 
seen in the dashed blue curve of the top frame of Fig. [1] 

This is to be contrasted with the case of /j, 7^ shown 
in the lower frames. Starting again from the left, the first 
arrow describes transitions from —62 +ei as before but 
its color has been changed to blue to indicate that now 
there is a restriction that only transitions with 17 > 'y+2fi 
are allowed. The next type of transition remains un- 
changed for < 7, as is the case shown. However, new 
transitions shown as the red arrow from +ei — *■ -1-62, can 
now occur, which were not possible for fi — 0. (The green 
shading represents the newly filled states due to finite /x 
which were formerly part of the empty conduction band.) 
Intraband transitions from +ei +ei are also possible 
and correspond to the delta function centered at 17 = 
which broadens to a Drude form when any kind of scat- 
tering is included. These transitions are depicted by a 
short red arrow with vanishing length for 17 — > 0. The 
next set from — ei +ei are restricted to 17 > 2/i while 
the last, —£2 ^ +£2, remain unaltered. The spectral 
weight lost between and 2fj, and between 7 and 7 -I- 2/t 
is balanced by the appearance of the two delta function 
contributions originating from the two new transitions 
+ei — *■ -1-62 and +ei +ei (Drude). The first set of 
these transitions leads to a delta function at 17 = 7 and 
this is traced to the perfect nesting of the bilayer dis- 
persion curves for ei (k) and £2 (k) which differ only by a 
constant displacement of 7. Should the electronic disper- 
sions stop being nested due to, for example, higher neigh- 
bor hoppings, or as we will see, due to the opening of the 
semiconducting gap, the delta function will broaden into 
a "band" , the width of which is related to the mismatch 
of the dispersions from perfect nesting. In particular, as 
we have already stated, the hopping from Al to B2 can 
be significant in size and will lead to changes in the band 
structure which are not cylindrically symmetric in the 
(kx, ky)-pla,iie^ . This fact complicates the mathematics 
and goes beyond the present work. However, the trigo- 
nal warping which arises can lead to subtle effects such 
as described by Mikitik and Sharlai^^ in relation to the 
Dirac points. 

If /i > 7, then further new transitions become possible. 
We now show an example for the case where /i = I.27 so 
that the -|-£2 band is now partially occupied. The possible 
transitions are shown in the bottom righthand frame of 
Fig. [T] For the first time all the terms in Eq. ([TO]) become 
activated. From the diagram —£2 —^ +ei transitions are 
restricted to 17 > 2^ + j (first term of second line in 
Eq. ([TO]) ) and the restriction on — £1 — > -|-£2 is 17 > 2(/i — 
7) -t- 7 = 2/t — 7 (second term of second line in Eq. (HH)). 
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FIG. 1: (Color online) Shown to the right. Upper frame: 
Frequency dependent conductivity a{Q,) of the bilayer nor- 
malized to twice that of a single graphene sheet ctq versus 
0/7 for three values of the chemical potential, as indicated, 
and A = 0. Middle frame: Dispersion of ei (solid line) and 
£2 (dashed) of the bilayer near the K point with the bands 
split by 7, showing the four types of transitions possible when 
^ = 0. Bottom frames: Transitions possible in the band struc- 
ture for finite chemical potential of = O.27 (lefthand side) 
and 1.27 (righthand side). See text for discussion. 



Both the — ei +ei and —€2 +£2 arise only for > 
2/i (first line of Eq. (|19p ) and new intraband transitions 
+€2 —>■ +£2 add the second term in a(/i). The b{fj,) is 
depleted through an additional negative contribution for 
fi > "f, because the +ei +£2 nested transitions are 
partially blocked by the filled states in +€2 near fc = 0. 
These facts manifest themselves in the top frame of Fig. [T] 
as a large peak centered at ft = and another at = 7 
from a(/i) and b{fi), respectively. All other transitions are 
completely suppressed below = 2^ = 2.47 except for a 
small contribution from — ei +62 clearly seen between 
2/x — 7 = 1.47 and 2.47. In addition, the added possible 
transitions above 2/i -I- 7 = 3.47 are seen to provide a 
small additional jump in a{Q) at this energy. 

In Fig. [31 we compare the results obtained on the 
basis of Eq. (|19p with the numerical evaluation of the 
complete formula Eq. fTB]) . For this evaluation a fi- 
nite value for the electron scattering rate l/rimp = '27] 
is needed. We took it to be a constant in energy and 
equal to 0.047. This is the simplest of models and is 
sufhcient for our purpose here. In general, the electron 
self-energy is a frequency-dependent complex function 
with non-vanishing real part. For impurity scattering, 
the detailed energy dependence of these functions will 
depend on the scattering potential and will be different 
for strong (unitary) and weak (Born) scattering as dis- 
cussed at length in a recent preprint^^. The energy de- 
pendence of the scattering rate can lead to interesting 
effects as discussed by Gusynin et ali^ for the case of 
the microwave conductivity where a cusp-like behavior 
is predicted to arise insead of the Lorentzian-like depen- 
dence of the usual Drude form. Here, it is sufficient to use 
a constant scattering rate in which case the real part of 
the self-energy, which is Kramers-Kronig related, is zero. 
Returning to the comparison in Fig. [Sj we note that for 
the analytic formula we broadened out the delta 
functions into Lorentzians using S{x) = {T/tt)/{x'^ + r^), 
where F = 1/Timp = O.O47 to match the full numerical 
work. The agreement between the numerical calculation 
(solid red curve) and the results of formula (dashed 
blue curve) is excellent. The differences arise solely be- 
cause no scattering was included in the dashed blue curve 
beyond broadening the two delta function contributions 
at = and = 7 for ease of comparison. We have 
verified that making 1 /rimp smaller brings the two curves 
closer, as it must. Note, that concurrent numerical work 
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to ours^^ shows the case of 2/x = 7 with varying impurity 
scattering. Their results agree with ours. Similar good 
agreement between our analytic results based on Eq. ([TO]) 
and the full numerical evaluation of Eq. (jl6p is seen in the 
lower frame of Fig. [3] for /i = I.27. Indeed, we see that 
broadening the delta functions of the analytical formula, 
as we have done, is essential to capture the result that 
the optical absorption below 2fj, — 7, in this case, is not 
zero but actually finite everywhere. This makes Eq. (flQl) 
a useful formula for experimentalists. 

The issue of optical spectral weight redistribution with 
changes in chemical potential can be addressed in a more 
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FIG. 2: (Color online) The evolution of the positive frequency 
spectral weight VF, found under each of the delta functions 
in Eq. (|19|) . as function of chemical potential. The case of an 
uncoupled graphene bilayer is given for comparison and the 
quantity c(/i) represents the spectral weight missing in the 
conductivity for < < 2/i relative to the /i = case. 



global fashion than we have done so far by introducing 
the partial optical sum: 



m) 



0+ o-Q 



(22) 



defined as the area under the conductivity up to en- 
ergy il. This is shown in Fig. U] for five values of the 
chemical potential: /1/7 = (solid black), 0.3 (long- 
dashed green), 0.5 (short-dashed-dotted red), 0.9 (long- 
dash-dotted pink) and 1.2 (short-dashed blue). In all 
cases by 17/7 = 3, the highest frequency shown, the in- 
tegrated spectral weight has returned to its /i = value 
(solid black curve) and the introduction of a finite charge 
carrier imbalance has not changed the partial optical sum 
up to that energy, although it has significantly changed 
its distribution in energy in the range < J7 < 87. Note 
in particular the sharp rise out of O = exhibited by 
all curves except for ^ = 0. This reflects the presence of 
the delta function contribution at 11 = which increases 
with increasing \x. The curves start to flatten when most 
of the spectral weight of the Drude contribution is inte- 
grated and consequently, this plateau occurs at about the 
same value as a(/i) plotted in Fig. [21 A second steep rise 
is also seen at 11 = 7 due to the second delta function 
in Eq. For the long-dashed green curve, we note 

the abrupt change in slope at H = 2/i = O.67 which has 
moved to 7 in the short-dash-dotted red curve, to I.87 
in the long-dash-dotted pink curve and to 2.47 in the 
short-dashed blue curve. These all reflect the 2/Lt cutoff. 
Note also in this last case, the small kink at 2/i — 7 which 
reflects the onset of the — ei +£2 optical transitions. 
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m/7=1.2 




4.0 



FIG. 3: (Color online) Comparison between the full numeri- 
cal evaluation of the conductivity [Eq. (I16|l ] and the analyti- 
cal formula of Eq. (|19p . For Eq. (|19p the delta functions are 
broadened into Lorentzians with a scattering rate to match 
the effective scattering rate in the numerical work. Two dif- 
ferent regimes of are shown: fi — O.27 (upper frame) and 
1.27 (lower frame). 



IV. RESULTS WITH ANISOTROPY GAP 

For the finite gap case, we begin with a discussion of 
the important energy scales involved. In the top frame of 
Fig. [5l we show the dispersion curves for our two bands: 
ei (solid red curves) and €2 (dashed blue curves). The 
first band shows a mexican hat structure. For positive en- 
ergies, there is a local maximum at zero momentum (p = 
Tivpk = with ti = A/2 = £^01 • There are also two min- 
imum at finite (f) = ±(A/2) ^(AM^VVlAM^T^) = 
fivpko with energy A7/(2-\/A2 -f 7^) = Egi. The low- 
est energy for which optical transitions are possible in 
the semiconducting case is 2Egi, when there is no charge 
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n/7 

FIG. 4: (Color online) Partial optical sum 1(0.) in units of 7 
versus n/7 for various values of chemical potential, as indi- 
cated in the figure. 

carrier imbalance which would introduce a finite /x. For 
understanding optical transitions, other energies are also 
significant and these are indicated in Fig. [S] At fc = 0, 
e2{k = 0) = ■\/yM-(A/2)^ = E02 and at the momen- 
tum of the minimum in band 1, the energy in band 2 



is £2(^0) = V(^'+7^ + 9AV/4)/(7' + A2) = Eg2. 
In addition, when a finite value of chemical potential 
is considered with fi > A7/ (2a/ -|- 7^) as shown, 
two other energies are important, namely the energy 
£2 for momenta at which /i crosses the ei dispersion 
curve. These are + A"^ + fi^ + 2L = Eb and if 
H < A/2, ^72 + A2 +^2_2L = Ea, where L = 
y/fi'^{A'^ + 7^) — (7A/2)2. These quantities determine 
the onset of various processes as we will describe below 
and are entered in Table HI Before doing so, it is impor- 
tant to understand how the electronic density of states 
N{e) varies with e since its value at the initial and the 
final value of the energy for a given optical transition pro- 
vides an important weighting factor for such processes. 

In the general case, one can obtain an analytic alge- 
braic expression for total double-spin N{e) in terms of the 
partial density of states provided by band 1 and band 2, 
A^i and N2, respectively, 

iV(e) = 7Vi(e) + 7V2(e), (23) 

where 

7Vi(e) = [7V+(e)-iVr(e)e(i?oi-e)]e(e-i?,i),(24) 
N2{e) = iV2-(e)e(e-So2), (25) 

with 



^;iVp i + (-i)"(7' + A2)/c" 



C = J 74 + (y + A2)[4e2 + A2 ± 2 V4e2(A2 + 72) - 72 A2] 



(26) 



The two terms for Ni simply refiect the two pieces of the 
energy dispersion for ei: the large k piece and the small 
k piece associated with the mexican hat. The density 
of states normalized by ^/h^vjp is shown in the bottom 
frame of Fig. [5] for several values of A. The long-dashed 
green curve is for the uncoupled graphene bilayer with 
7 = A = and is included for comparison. We note that 
beyond e/7 ^ 1.2 for the parameters shown here, the var- 
ious curves come close together on the scale of the figure. 
The inset show the density of states associated with each 
band separately for the case of A — I.57. The ei dis- 
persion gives rise to the square root singularity, followed 
by a shoulder, that is seen in the total density of states. 
Indeed the singularity in both the solid and short-dashed 
curves of the main frame at e = Egi and the shoulder at 
A/2, have their origins in the mexican hat structure of 
the dispersion curve. The shoulder comes from the top 
of the hat and is a van Hove singularity associated with 
the dispersion flattening at this point. The square root 



singularity can be derived for e near the hat minimum 
e ~ Egi and we find, as othersS^ have done, 

kn 2m* 

^W = t^a/ E^- 27 

' AttH^ e- Egi 



The factor kg = (A/2?iwf)\/(A2 + 2j^)/{A^ + 72), 
which is the momentum associated with the rim of the 
mexican hat, is related to the degeneracy of energies 
around the circular minimum of the mexican hat (fixed 
magnitude, varying angle of momentum) while the square 
root is associated with the one-dimensional variation in 
energy up and down the rim with varying magnitude of 
k (fixed angle), making this equivalent to what is ex- 
pected for a one-dimensional density of states. The ef- 

r_^2 I 2i3/2 

fective mass is m* = 2a'd^ a^+2j'^ with energy e — 
Egi + 2^(fc — ko)^. For comparing with Eq. it is 
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TABLE I: Energies involved in optical transitions (see top frame of Fig. [5]) for the case 7 — 0.4 eV, A — 0.6 eV and two values 
of chemical potential /j, = 0.25 and 0.35 eV. 





Formula 


value in oSl 


value in eV 






for /i = 0.25 eV 


for /i = 0.35 eV 




7A/(2V7' + A2) 


0.166 


0.166 


Eg2 


V(A4 + 7" + 9A27V4)/(72 + A2) 


0.74 


0.74 


Eol 


A/2 


0.3 


0.3 


E02 


Vt^ + (A/2)2 


0.5 


0.5 


Ea 


-^72 + A2 + - 2i 


0.56 




Eb 


^^2 + A2 + ^2 + 2i 


0.92 


1.042 




where L = a/m2(^2 j_ ^2) _ (7A/2)2 







convenient to rewrite Eq. (I27p in the form: 



NU) 



'7A^A2 +72 



E, 



91 



(28) 



This singularity plays an important role in optics because 
it leads to peaks in a{Q) when the energy of either or both 
the initial and final states involved falls at Egi. This fact, 
plus the energy scales identified in Fig. [5] (top frame) 
allows one to understand the qualitative features of our 
numerical results for the conductivity. 

Referring to Fig. [6l we start with the upper frame, 
where A = O.27. Three values of the chemical poten- 
tial fj, are considered. The dashed blue curve, which is 
for reference, has = 0, the dash-dotted black is for 
fi — O.I7 and the solid red is for fj, = O.27. We be- 
gin with the dashed blue curve which is the only case 
which exhibits a true semiconducting gap with no ab- 
sorption up to photon energy Q = 2Egi = O.I967. We 
first note the existence of a small tail below this energy 
and also note that the main rise is smeared. This is be- 
cause we have used a finite impurity scattering rate of 



imp 



= 0.047 in all of our numerical work. It is clear 



that most of the states lost below the gap {2Egi) are to 
be found piled up in a region of order a few Egi above 
it. Later we will examine in detail this optical spectral 
weight redistribution. This part of the curve for a{Q) is 
due to the — ei +ei transitions and their onset involves 
the singular density of states at Egi for both initial and 
final states. The second peak in the dashed blue curve is 
traced to transitions —62 +ei and — ei — > +€2 (degen- 
erate in energies) and these involve a single singular N{e) 
either in the initial or final state but not for both, and 
the onset is at E02+E01 = + (A/2)2 + A/2 ~ I.I7. 
It is also clear from Fig. [5] (top frame) that this onset re- 
mains even for finite fi because the —ei—f +62 transition 
is never blocked by a finite fi value provided /z falls below 
the minimum in £2 which is the case considered here. 

Next we consider the dash-dotted black curve for 
fi = O.I7. In this case, the — ei ~f +ei transitions 
are no longer possible for energies less than or equal to 
2/i = 0.27 because these states are occupied and cannot 
be used as final states in optical transitions. The energy 
2/i = 0.27 is very close to the onset for the dashed blue 



curve with ^ = which is at O.I967 and, as a result, both 
dashed and dash-dotted curves are close to each other in 
this energy region. Nevertheless it is important to un- 
derstand that while the dashed blue curve is in principle 
encoded with the information on the details of the mex- 
ican hat topology, the dash-dotted black curve is less so 
as it involves, in addition, a sharp cut off at O.27, which 
obscures some of these details. Finally we note the trans- 
fer of optical spectral weight to a Drude peak centered 
at Q = due to +£i — > +ei intraband transitions. This 
peak does not exist at zero temperature for A 7^ and 
/i = 0. 

The next interesting feature in the dash-dotted black 
curve is the peak slightly above fl — 0.9j. This peak can 
be traced to new +ei — > +£2 optical transitions made 
possible for finite /i. We have already seen that, for 
A = 0, these would fall at O = 7 and provide, in the clean 
limit, a delta function contribution which can be traced 
to the perfect nesting of the dispersion curves. But for 
finite A, perfect nesting no longer occurs and these tran- 
sitions broaden into a "band" as well as shift in energy 
as we now describe. Returning to Fig. [5] (top frame), we 
note that because +ei states are now occupied in a re- 
gion about the mexican hat minimum, transitions from 
+ei — > +£2 become possible for photon energies between 
Ea—^^ and Es—fJ^- The intermediate energy transition at 
Eg2 — Egi involves the density of states at the bottom of 
the mexican hat and this is expected to lead to a peak in 
a{il) at this energy. For the parameters of the model, the 
new "band" ranges in energy between 0.9157 to 0.9347, 
with a peak at O.937. In Fig. [51 an additional broadening 
of this narrow "band" is included because we have used 
a finite l/rimp- Nevertheless, overall, the new "band" 
does not broaden much as a result of scattering and the 
changed topology of the energy dispersion curves ei(k) 
and e2(k). While the dispersion curves +ei and -1-62 are 
no longer simply displaced by a constant amount with re- 
spect to each other and therefore are no longer perfectly 
nested, the effect is not large. 

As we have already discussed, the — ei — > +62 transi- 
tions remain unaffected by a finite but small value of fi 
and the onset for these transitions remains unshifted in 
energy at Eqi + E02 ~ 1.1. However, the optical weight 
at the onset is depleted because the —62 +ei transi- 
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FIG. 5: (Color online) Top frame: Band structure around the 
Fermi level in the presence of the asymmetry gap A for real- 
istic values of A = 0.6 eV and 7 — 0.4 eV. A finite chemical 
potential of fi = 0.25 eV is shown. Various important energies 
are indicated in the figure and displayed in Table |I1 Bottom 
frame: The density of states N{e), in units of j/h^vp, for 
the bilayer for several values of A as indicated in the figure, 
in comparison with the case of 7 = A = corresponding to 
the uncoupled graphene bilayer. The inset shows the partial 
density of states A'^i (e) and N2 (e) for the two separate bands 
ei and £2, respectively, for the case of A/7 = 1.5. 



tions do know about ^. For /x = O.I7, these transitions 
are blocked in the energy range Ea + /i = 1.1 167 to 
Eb + M = 1.1347. Both these energies fall very close to 
the peak energy I.I7 and show up in the figure simply as 
a slight depletion of the large broad peak above = I.I7 
in the dashed blue curve. 

Similar arguments explain the main features seen in 
the solid red curve for /z = O.27. In this case, the cutoff 
imposed on the — ei — s- +ei transitions is at $1 = O.47 
{2fi) and the amount of spectral weight transferred to 



FIG. 6: (Color online) Optical conductivity in the presence of 
an asymmetry gap A. The upper frame shows the results for 
varying /i with fixed A = O.27. The bottom frame shows the 
result for realistic values of the parameters, i.e., A = 0.6 eV 
and 7 — 0.4 eV. Here two curves are shown, one for /x = 0.25 
eV where the chemical potential lies above the gap in the 
band structure but below the hat maximum, and fi — 0.35 
eV, where it lies above the hat maximum. 



the Drude centered at fl = has greatly increased as 
compared with the dash-dotted black curve. The peak 
from the new +ei —>■ +62 transitions has broadened ex- 
tending from 0.97, set by Eq2 — Eqi (as /i is now above 
the hat maximum) to 0.9987, set by Eb — fi- Finally, the 
upper cutoff on the —62 — s- +ei transitions has moved to 
Eb + M = 1-4 eV, where a sharp rise in conductivity is 
seen. 

In a field effect device with a semiconducting gap, the 
value of A cannot be set independently from the value 
of the chemical potential which sets the occupation of 
the +ei band and Egi, the energy of the minimum of the 
mexican hat. Realistic values based on numerical work 
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presented in Ref;^ are 7 = 0.4, A = 0.6, and n = 0.35. 
Results for the conductivity (7{fl)/(7Q versus are shown 
as the dashed red curve in the bottom frame of Fig. [6l 
We found it illuminating also to present for comparison 
the case oi fi — 0.25 represented by the solid blue curve. 
The band structure involved is shown in the top frame of 
Fig. [51 where various energies of importance for optical 
transitions are identified. Actual numbers for the i?i's are 
found in Table HI All the important qualitative features 
seen in the conductivity curves can be understood on the 
basis of these energy scales. We begin with the solid 
blue curve. If the chemical potential was zero the — ei — > 
+€i transitions would start at 2Egi = 0.333 eV but for 
^ = 0.25 these are cutoff at 0.5 eV where we see a sharp 
rise in conductivity. On the other hand, the +ei — > +62 
transitions, made possible through the finite occupation 
of the +ei states at finite chemical potential /i, start at 
Ea — fJ, = 0.31 eV where we see the first rise in a{^l) after 
the Drude. These transitions extend to Eb — IJ^ = 0.67 
eV where the conductivity shows a steep drop. Thus, the 
nested transitions of the A = bilayer graphene case are 
now greatly broadened and their onset has moved below 
the onset for the — ei +ei transitions. The small peak 
at Eqi + E02 = 0.8 eV is traced to the —62 +ei and 
— ei +£2 transitions which, as we have discussed, are 
depleted by charging but their onset remains unshifted. 
Finally the rise ati?B+/j,= 1.17 eV coincides with the 
upper cutoff on the —&2 ^ +ei transitions. For ji ~ 
0.35 eV (dashed red curve) the +ei +€2 transitions 
start at i?02 — £^01 = 0.2 eV and cause the first rise in 
conductivity which follows the Drude peak about il = 
0. They extend to up Eb — h = 0.69 eV where the 
conductivity shows a rapid decrease. We stress that the 
onset of these transitions is set by the difference {Eq2 — 
Eqi) (bottom of second band minus top of mexican hat) 
and not directly by the chemical potential which falls 
above Eqi. The peak at 0.8 eV is further depleted as 
compared with the solid blue curve but is unshifted. The 
upper cutoff on the —62 — > +ei transitions is now at 
Eb+ 1.39 eV. 

In Fig. [71 we show results for the partial optical sum 
/(II) in units of 7 versus fl/^ which describes the opti- 
cal weight redistribution brought about by the opening 
of the anisotropy gap and finite chemical potential. Four 
cases are compared in the main frame for /i = 0. The 
dash-dotted green curve is for A = 0, long dashed black 
for A/7 = 0.5, solid red for A/7 = 1.0, and short dashed 
blue for A/7 = 1.5. We see the almost complete deple- 
tion of optical spectral weight below the gap 2Egi set 
by the minimum energy in the mexican hat dispersion 
curves. The small tails below this cutoff are due to the 
impurity smearing caused by our use of a finite 1 / Timp in 
the numerical work. In all cases, the lost spectral weight 
is recovered in the region above it. The overshoot be- 
yond the reference curve (dash-dotted green) is largest 
for the largest A considered and in this case /(f^) is still 
slightly above the A = reference even at = 3 eV. It 
is clear, however, that there is a close balance between 
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FIG. 7: (Color online) The partial optical sum /(fi) in units of 
7 versus $7/7 for various values of A with ^ = 0. Inset: li^C) 
for the more realistic parameters used for (j(f2) in Fig.[6]D, i.e., 
7 = 0.4 eV, A = 0.6 eV, and = 0.25 eV (dashed red curve) 
and 0.35 eV (solid blue curve). 



weight lost and gained in this interval. In the inset of 
Fig. [71 we show /(r2) in units of 7 versus 5^/7 for more 
realistic values of 7, A, and /i associated with the con- 
ductivity curves shown in the bottom frame of Fig. [51 
For [I finite, there is a narrow Drude peak in (T(f2) with 
width set by l/rimp, which results in a fast rise of /(fi) 
out of r2 = 0. This is followed by further rises with small 
kinks reflecting the sharp rises and drops in the associ- 
ated conductivity curves. 



V. SUMMARY AND DISCUSSION 

We have derived an analytic algebraic formula for the 
zero temperature, clean limit optical conductivity of a 
graphene bilayer valid for any value of charge imbalance 
characterized by a general value of the chemical potential 
/i. In the limit of /z = our formula reduces properly to 
that of Abergel and Fal'ko^. A finite /x introduces sev- 
eral important modifications. First, some of the low en- 
ergy — ei — > -|-ei transitions are no longer allowed because 
the charge imbalance leads to finite occupation of part of 
the -|-ei band. Furthermore, new — s- +ei intraband 
transitions are now possible and this leads to a finite 
Drude peak at i7 = 0. In addition, new interband tran- 
sitions between +ei — > ^t2 are allowed and these lead 
to a second delta function contribution at 17 = 7. These 
transitions reflect the perfect nesting between bands 1 
and 2 which are simply displaced in energy by a constant 
amount 7 in our model band structure when there is no 
anisotropy gap. For more complicated models in which 
the nesting ceases to be perfect, the peak at f2 = 7 would 
broaden into an absorption band whose width and struc- 
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ture in energy reflects the mismatch in topology between 
bands 1 and 2. In addition, there is a broadening brought 
about by the elastic scattering rate l/rimp- Immediately 
above the = 7 peak there is a second depletion region 
which is due to the blocking of the — £2 — *■ +ei transi- 
tions between 7 and 7 + 2/i. The presence of a finite fi 
does not affect the — ei — > +€2 transitions which also fall 
in the same energy interval and so only partial depletion 
is involved. It is found that the optical spectral weight 
lost in the two depletion regions is completely compen- 
sated for by the spectral weight which resides in the delta 
function. We have just described the case for /i < 7. For 
fi > J, new intraband transitions associated with the 
+62 +62 transitions become possible and the transi- 
tions which utilize band 2 as the final state also become 
blocked. 

Note that all the main features of the conductivity 
curves found here can be traced directly to the un- 
derlying band structure of bilayer graphene. Correla- 
tion effects, which have not been treated here, lead to 
self-energy (E) corrections. The real part of S renor- 
malizes the single particle energies and the imaginary 
part introduces damping. Angle-resolved photoemission 
spectroscopy (ARPES) provides a direct measure of the 
dressed dispersion curves and of their many body broad- 
ening. The ARPES data^"^- confirms the general shape of 
the tight-binding dispersion curves used here with some 
smaller modification and by implication, we do not ex- 
pect large changes in the frequency dependence of the 
conductivity described here. Of course, the scattering 
rate due to many body interactions will, in general, be 
frequency dependent while here we have treated it as a 
constant. In bilayer graphene, even impurity scattering 
involves an energy-dependent scattering rate and many 
examples of how this changes the shape of the optical 
conductivity can be found in the recent work by Nilsson 
et al.-^ Finally, we note that vertex corrections have not 
been included here as this goes beyond the present work. 

When a semiconducting gap is introduced, the band 
structure becomes modified. In particular, band 1 ac- 
quires a mexican hat structure with the top of the hat 
at energy Eqi = A/2 and the minimmn on the rim of 

the hat is at Egi = {A-f /2) / y/^^Tl^ ■ This can be 
achieved in a graphene bilayer when donor atoms are 



seeded on its top surface and the whole is placed in a 
field effect configuration. In such junctions, finite A also 
implies finite ^ > Egi and there exists a finite charge 
imbalance. Only for the case of finite A and ^ = is 
there a real gap in the system, i.e. no absorption below 
2Egi: the — ei —>■ +ei transitions become gapped and 
the lost optical spectral weight is found to accumulate 
in the energy region just above the gap on an energy 
scale of order A. In any practical case, however, a fi- 
nite value of chemical potential accompanies a finite A 
and there is a Drude centered at = and a cutoff 
on the — ei +ei transition of 2^. The onset of the 
— ei +62 (and —£2 — *■ +£1) moves from 7 in the case 
of A = to £^01 + £^02 = VW2F+7^ + (A/2). The 
—62 +ei transitions become gapped in the interval 
Ea + /i to Eb + IJ., unless /i > A/2, then the lower limit 
changes to Eqi + Eq2. In addition, the transitions from 
+ei +£2, which are possible for finite fi and which in 
pure unbiased bilayer graphene provide a delta function 
at 7, are broadened into a band from Ea — /i to Eb — /i, 
the width of which depends on the value of chemical po- 
tential. This band of absorption starts at energies below 
7, with onset at the energy of the top of the mexican hat 
if n is greater than this energy. While these modifications 
in the possible transitions due to finite A and fj, can lead 
to complicated spectral weight shifts in a{n) versus f2, 
all changes can be understood qualitatively from a knowl- 
edge of the band structure involved. Biased graphene bi- 
layers offer a rich pattern of variation of a^fl) versus O as 
the size of the anisotropy gap is varied through changes 
in the voltage of the field effect device. This system is 
not only important because of possible practical appli- 
cations, but it is also the only system known for which 
the value of the semiconducting gap can be tuned by the 
application of an external voltage. 
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